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Abstract. We briefly summarize the method of simulating Sgr A* polarized sub-mm spectra 
from the accretion flow and fitting the observed spectrum. The dynamical flow model is based on 
three-dimensional general relativistic magneto hydrodynamic simulations. Fully self-consistent 
radiative transfer of polarized cyclo-synchrotron emission is performed. We compile a mean sub- 
mm spectrum of Sgr A* and fit it with the mean simulated spectra. We estimate the ranges 
of inclination angle 8 — 42° — 75°, mass accretion rate M = (1.4 — 7.0) x lO~ 8 M0year _1 , and 
electron temperature T e = (3 — 4) x 10 10 K at 6M. We discuss multiple caveats in dynamical 
modeling, which must be resolved to make further progress. 
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1. Estimating Sgr A* BH and accretion flow parameters 

Sagittarius A* (Sgr A*) is a supermassive black hole (BH) in the center of the Milky 



Way galaxy Its accretion flow is underluminous and radiatively inefficient (|Yuan et al 



2003; Sharma et al.ll2007l ) for the radial profiles of electron temperature T e consistent with 
observations. Thus the gas infall can be modeled with three-dimensional (3D) general rel- 
ativistic (GR) magneto hydrodynamic (MHD) simulations without c ooling. Such simula- 



tions and the associated rad i ative transfer were rec ently performed by Moscibrodzk a et al 
(l2009l) : lDexter et al.l (|2010l) : Iffilburn et all (|2010h and by our group ([Shcherbakov et al 
2012h ; while some groups re sorted to (semi-)analytic accretion flow models (jHuang et al 
2009; Broderick et al.ll2011 ). All researchers aimed at fitting combinations of radio, (po- 
larized) sub-mm, IR, and X-ray spectra and emitting region size. 

Our group adopts the following approach in fitting Sgr A* observations and estimating 
the accretion flow and the BH parameters. We compile the mean Sgr A* polarized sub- 
mm spectrum from observations in the past ~ 10 years. We perform a set of 3D GRMHD 
simulations without cooling of accretion flow onto the BH for several spin values. We ap- 
ply the original GR polarized radiative transfer code, which follows the prescription in 
Shcherbakov fc Huand (l201lh and employs Faraday rotation and conversion coefficients 
from Ishcherbakovl (|2008l ) . The radiative transfer is fully self-consistent. The mean simu- 
lated spectrum is found as a mean of spectra at different simulation times. We explore the 
parameter space of spin a*, inclination angle 9, mass accretion rate M , and the ratio of 
proton to electron temperature T p /T e at radius 6M. For each spin we find a model, which 
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provides the best agreement between the mean simulated and the mean observed spectra. 
We test our numerical radiative transfer code and the results for robustness. Our best 
bet model has a* = 0.5, 6 = 75°, M = 4.6 x lO^Mgvear" 1 , T e = 3.1 x 10 10 K at 6M, 
and spin position angle PA = 115°. Model parameters are in the ranges 6 = 42° — 75°, 
M = (1.4 - 7.0) x 10- 8 M Q year-\ and T e = (3 - 4) x 10 10 K at 6M across all spins. 
The sub-mm circular polarization is mainly produced by Faraday conversion as modified 
by F araday rotation. The emission region size at 230 GHz is consistent with the VLBI 
size (Doeleman et al. 2008T) . The details of the modeling can be found in our main paper 
([Shcherbakov et alJ^oM T 



2. Caveats 

Present-day dynamical modeling of tenuous plasma falling onto the BH inevitably 
relies on heavy approximations. Our numerical resolution in MHD simulations might 
not be s ufficient to ach i eve co nvergence in turbulence dissipation rate or magnetic field 
structure lHawlev et al. ( 2011 ). The dynamic range and duration of the simulation are not 
large enough to obtain reliable radial profiles of quantities and the long-term behavior 
of the accretion flow. MHD approximation itself might not hold for collisionless plasma. 
Non-trivi al plasma effects include n on-thermal distrib ution of particles inferred for both 
the flared Dodds-Eden et al. ( 20091 ) and the quiescent Yuan et al. ( 20031 ) states. Another 
important effect is thermal heat con duction. It was found to dramatically alter the radial 
densi ty profile and other quantities ([Johnson fc Quataertll2007l : IShcherbakov &: Baganofl 
2010). In view of the aforementioned difficulties the results of our modeling should be 
viewed as estimates with unknown systematics. We have developed the framework of 
testing various dynamical models of Sgr A* and other low-luminosity active galactic 
nuclei. Whenever the more self-consistent dynamical models become available, it will be 
possible to promptly apply them to fit Sgr A* accretion flow. 
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